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Abstract. We report the first simultaneous measurement 
of the broad band X-ray (0.3-150 keV) spectrum of the 
neutron star x-ray binary 4U 0614+091. Our data con- 
firm the presence of a hard x-ray tail that can be mod- 
eled as thermal Comptonization of low-energy photons on 
electrons having a very high temperature, greater than 
220 keV, or as a non-thermal powerlaw. We detected 
a spectral feature that can be interpreted as reprocess- 
ing, via Compton reflection, of the direct emission by an 
optically-thick disk and found a correlation between the 
photon index of the power-law tail and the fraction of radi- 
ation reflected which is similar to the correlation found for 
black hole candidate x-ray binaries and Seyfert galaxies. 



Key words: accretion, accretion disks — stars: individual 
(4U 0614+091) — stars: neutron — X-rays: stars 

1. Introduction 



The x-ray source 4U 0614+091 was discovered in 
the Uhuru survey and was localized by Copernicus 
(Willmore et al. 1974). The source is an x-ray burster 



( Swank et al. 1978] ) and, thus, contains a neutron star. 
The luminosity of the x-ray bursts constrains the distance 



to be within 3 kpc (Brandt et al. 1992). Classification of 
4U 0614+ 091 as an atoll source was suggested by Si ngh fc 
Apparao (1994) and confirmed by Mendez et al. (1997). 
New study of the source has been stimulated by the discov- 
ery of high frequency quasi-periodic os cillations (QPOs) 
in the persistent emission of this source ( Ford et al. 1997 ) 
and of a hard x- ray tail extending to at least 100 keV 
( |Ford et al. 1996]) . 

Here, we report on observations of 4U 0614+091 ob- 
tained with the narrow fi eld instruments (N FIs) on board 
the BeppoSAX satellite ( |Boella et al. 1997|) . We have ob- 
tained the first simultaneous measurement of the broad 



band x-ray spectrum of 4U 0614+091, extending from soft 
x-rays to hard x-rays (0.3-150 keV) and are able to place 
strong constraints on the spectrum of the hard x-ray 
emission. We find that the continuum hard spectrum of 
4U 0614+091 is well described by a powerlaw or thermal 
Comptonization with a reflection component. We find that 
the thermal Comptonization temperature and the lower 
bound on any exponential cutoff of the powerlaw model 
are both greater than 200 keV. 

We also studied the relation between the reprocessed 
component and the intrinsic component. A correlation be- 
tween the strength of reflection and t he intrinsic spec- 
tral slope has been found for GX 339-4 ( Ueda et al. 1994 ) 
and recently for Seyfert AGNs, 2 black holes candidate 
(BH C) X-r ay binaries, and 2 X-ray bursters (Zdziarski 
et al. |1999| ). We find a similar relation for 4U 0614+091. 
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2. Observations and Analysis 



We performed a joint BeppoSAX/RXTE observation 
of 4U 0614+091 on 19-20 October 1998 for a total 
of 42.6 ks of on-source observing time in BeppoSAX 
( Piraino et al. 199S| ). Here, we use data from the four Bep- 
poSAX Narrow Field Instruments (NFIs) in overlapping 
energy bands selected to give good signal to noise for 
this source: the Low Energy Concentrator Spectrometer 
(LEGS) for 0.1-4 keV, the Medium Energy Concentrator 
Spectrometer (MEGS) for 1.8-10 keV, the High Pressure 
Gas Scintillation Proportional Counter (HPGSPC) for 4- 
27 keV, and the Phoswich Detection System (PDS) for 
12-200 keV. LEGS and MEGS data were extracted in cir- 
cular regions centered on the source position using radii 
of 8' and 4' respectively, containing 95% of the source 
flux. The spectra have been rebinned to have at least 30 
counts per channel, and the HPGSPC and PDS spectra 
were grouped using a logarithmic grid. A normalization 
factor has been included to account for the mismatch in 
the BeppoSAX instruments' absolute flux calibration. The 



S. Piraino et al.: BeppoSAX Observations of 4U 0614+091 



fit values of relative normalization are in good agreement 
with values typically observed (Cusumano ct al. 1998). 

For the total spectrum presented in Figures |l| and g, 
we summed all of the 42.6 ks of data obtained in Octo- 
ber 1998. To study the correlation of spectral slope with 
reflection, presented in Figure |3[ we divided this obser- 
vations into two intervals, based on luminosity, and also 
used data from an observation performed on 13-14 Oc- 
tober 1996 available from the BeppoSAX archive, again 
divided into two intervals. 



3. Spectral Results 

We fit the 0.1-200 keV BeppoSAX spectrum of 
4U 0614-1-091 with several single component con- 
tinuum models: powerlaw, powerlaw with an expo- 
nential high energy cutoff, three different Comp- 
models 



tonization models (Sunyaev & Titarchuk 1980; 

[Titarchuk 1994 ; Poutanen & Svensson 1996), an ex- 
ponentially cutoff powerlaw and a Comptonization model 
( Poutanen fc Svensson 1996|) each with added reflection 
(Magdziarz & Zdziarski 1995), thermal bremsstrahlung, 
several disk models, and each of these single component 
models with an added single temperature simple black- 
body component. Each model also included absorption 
and a Gaussian line near 0.7 keV ( phristian et al. 1994 ; 
White et al. 1997t [Schulz 1999| ). Only the powerlaw (with 
or without an exponential cutoff) and Comptonization 
models gave close to acceptable fits {xl —1.33) and 
in each case the fit was improved by addition of the 
blackbody component. In Table 1 we report the best fit 
parameters for the powerlaw model and the best fitting 
Comptonization model (Poutanen & Svensson 1996) 
both with and without reflection. The fit residuals of 
the models without reflection show a feature between 10 
and 60 keV that disappears when reflection is added, see 
Figure w. The two models with reflection provided the 
best fit with a significant improvement in xt (^^^ Table 
1) relative to any other model. 

The powerlaw with reflection model is the 
PEXRAV/PEXRIV model in XSPEC, which is an 
exponentially cut-off powerlaw spectrum with reflection 



(Magdziarz & Zdziarski 1995) from a disk. We assumed 
solar abundances and fixed at 60° the disk inclination 
angle. Adding an exponential cutoff does not improve 
the fit, either with or without reflection, and the cutoff 
energy is well above 200 keV. The Comptonization model 
with reflection is COMPPS, kindly supplied by Juri 
Poutanen. This code calculates radiative transfer and 
Comptonization in a two-phase disk-corona geometry 
( Poutanen fc Svensson 19961). Reflection is calc ulated via 
the same method (Magdziarz & Zdziarski 1995) as for the 
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PEXRAV model. We used the model to calculate thermal 
Comptonization and used an approximate treatment of 
radiative transfer using the escape probability for a sphere 
for the results presented in Table 1. The temperature 



1 10 IOC- 

Channel energy (keV) 

Fig. 1. Residuals for the models described in Table 1. For 
display, the energy bins are equal to the la energy reso- 
lution. 



of the input soft photon distribution is driven towards 
unexpectedly low values when left as a free parameter. 
We tried fixing this temperature to that of the lowest 
temperature soft photon source present in our spectral fit 
(the 0.27 keV blackbody), but were unable to obtain good 
fits. For both reflection models, the ionization parameter 
is consistent with zero, but is not well constrained. For 
the results presented in Table 1 and Figures 2 and 3, 
we used the PEXRAV model, which assumes a neutral 
reflecting medium, with no exponential cutoff. Figure 2 
shows the BeppoSAX data and the various components 
of the fit using the PEXRAV model. 

As the presence of an absorption edge is an inher- 
ent feature of reflection, before accepting the reflection 
as component of the model, we verifled the presence in 
the spectrum of an absorption edge. Adding this feature 
to the powerlaw or Comptonization models, we found an 
improvement in the x^ and the same edge energy and op- 
tical depth (i?cdgc — 8.5 keV, r = 0.1) found by Singh & 
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Table 1. Parameters for fits of the spectrum of 
4U 0614+091 



1 10 100 

Channel energy (keV) 

Fig. 2. LEGS, MEGS, HPGSPG and PDS unfolded av- 
eraged spectra, for the October 19-20, 1998 observation 
together with a model consisting of a powerlaw, reflec- 
tion, a blackbody, and a low energy Gaussian line all with 
absorption. The total model fit is shown as a solid line, 
the powerlaw as dotted line, the reflection component as 
a dot-dot-dashed line, the blackbody component as a dot- 
dashed line, and the Gaussian line as a dashed line. 



Apparao (1994) in a low state observation of 4U 0614+091 
with EXOSAT. 

Adding a Gaussian emission line to model iron K^ flu- 
orescence did not improve the fit significantly and we place 
upper limits of 20-60 eV on the equivalent width for lines 
in the 6.4-8.1 keV range. While iron Kq fluorescence is 
commonly associated with reflection in AGN, the iron line 
emission found in Galactic BHG X-ray binaries in the low 
(hard) state is much less than that found in AGN (e.g. 



Zdziarski et al. 1998). This may be due to differences in the 



ionization state or geometry of the disk, Doppler broaden- 
ing and relativistic changes in the line profile which make 
t he line difficult to detect, or resonant Auger destruction 
( Ross, Fabian fc Brandt 1996 ). All of these effects may 
also hold for neutron star X-ray binaries. In addition, the 
equivalent width of iron K^ for 4U 0614+091 is likely to 
be low since the equivalent width of iron K^ decreases 
approximately linearly with photon index F as result of 
the decrease in the number of incident photons with en- 
ergy larger than the absorption edge energy for higher F 
( [George fc Fabian 199l| Vrtilek et al. |1993| ) and F is larger 
for 4U 0614+091 than for BHG x-ray binaries in the low 
(hard) state. 

We found a blackbody component with temperature, 
fcTbbj near 1.5 keV. As the Gomptonization optical depth 
is small, the detection of the blackbody is consistent with 
an origin within the Gomptonization region. We note 
that the blackbody flux found using GOMPPS is twice 
that found using the powerlaw. Thus, it is unlikely that 
the blackbody component is a spurious feature due to 
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Comp 


PL 


Comp 








with reflection 


Nh 


3.3+0.1 


3.0+0.1 


3.7+0.2 


3.3+0.3 


E, 


0.71+0.03 


0.70+0.03 


0.66+0.04 


0.65+0.05 


a 


0.07+0.02 


0.07+0.02 


0.09+0.02 


0.11+0.02 


EW 


88+33 


90+20 


140+64 


140+60 


r 


2.33+0.02 




2.44+0.03 




kT, 




<0.04 




<0.03 


kn 




160+35 




26011° 


T 




0.26+0.1 




0.10+0.03 


Refl 






0.90+0.25 


0.45+0.25 


kT^y, 


1.45+0.08 


1.28+0.15 


1.47+0.06 


1.39+0.06 


Fbb 


0.51+0.10 


0.75+0.15 


0.82+0.12 


1.5+0.3 


F 


16.73 


17.83 


16.86 


16.94 


Fun 


40.33 


38.39 


47.93 


34.17 


xi{y) 


1.21 (561) 


1.17(559) 


1.09 (558) 


1.08 (555) 



Note — Given in the table are the absorption column den- 
sity (A'^h) [in units of 10^^ cm"^], the centroid (Sg) [keV], 
width (cr) [keV], and equivalent width {EW) [keV] of the 
low energy emission line, the photon index (F) of the pow- 
erlaw, and the temperature of the input seed photon distri- 
bution (fcTs) [keV], the temperature (fcTo) [keV] and optical 
depth (r) of the Comptonizing electron cloud, relative reflec- 
tion fraction {R), the temperature (fcTbb) [keV] and flux (-Fbb) 
[10~^° ergcm"^ s~^] of the blackbody, the absorbed {F) and 
unabsorbed (Fun) flux [10"^° ergcm"^ s"^ for 0.1-200 keV], 
and the reduced chi square, x^, and degrees of freedom (j/). 
All quoted errors represent 90% confldence level for a single 
parameter. 



approximation of the Gomptonization spectrum with a 
powerlaw. Blackbody components have been found from 
4U 0614+091 in the past, with kT ranging from 0.3 keV 
( iGhristian et al. 1994|; ^chulz 19991) to 1 .7 keV (Singh & 



Apparao 1994; Barret fc Grindlay 1995). Adding a sec- 



ond blackbody component to our model, we found that 
the second blackbody improves the fit slightly and has 
kT = 0.27 ± 0.03 keV. Thus, two blackbody components 
may be simultaneously present in the spectrum of this 
source. Addition of the second blackbody does not sig- 
nificantly affect the other fit parameters. We also tried 
using the blackbody accretion disk models of Shakura 
fc Sun ayev ( |1973| ), Stella & Rosner ( |1984| ), and Mitsuda 
(1984). These models also gave good fits, but not better 
than the fits obtained with the simple blackbody. 

Using the PEXRAV model in a total of four intervals 
over the two observations, we found that the amount of 
refiection and also the blackbody flux are correlated with 
the photon index (Fig. 13) . The values found for the reflec- 
tion coefficient are large. This may indicate that the pho- 
ton flux irradiating the disk is larger than the observed 
flux as expected for relativistic electrons in an anisotropic 
refiection geometry (GhiscUini et al. 1991) or may simply 
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Fig. 3. Spectral parameters correlations, (a) Magnitude 
of reflection versus photon index, (b) Blackbody flux 
(erg cm~^ s~^) versus photon index. 



indicate that the inclination is less than 60°. The unab- 
sorbed luminosity in the 0.1-200 keV band varied over the 
range 5.1 — 29.1 x 10^^(3^-^)^ ergs"^ in the four intervals. 

4. Discussion 

Before the BeppoSAX era, there were no x-ray burster 
observations with simultaneous energy coverage from 0.1- 
200 keV. This broad coverage is very important to ac- 
curately measure the x-ray spectrum. Using BeppoSAX, 
we find a spectrum for the x-ray burster 4U 0614-1-091 
which can be described either with a powerlaw with no 
detectable cutoff below 200 keV or with a thermal Comp- 
tonization spectrum with an electron temperature in ex- 
cess of 220 keV. As equally good fits are obtained in the 
two cases, we cannot determine whether the emission is 
due to a thermal or non-thermal process. If the emission 
is due to Comptonization, with either a thermal or non- 
thermal electron energy distribution, then the input soft 
photon distribution must peak at an energy well below 
that of any observed soft photon source; in particular, well 
below the temperature of the cooler (0.27 keV) blackbody 
component. Production of these soft photons must be ad- 
dressed in any successful application of a Comptonization 
model to 4U0614-h091. 

The only robust evidence for the presence of a black 
hole in an X-ray binary is a measurement of the mass 
of the compact object indicating a value greater than 
3M0. While neutron stars have unambiguous x-ray ob- 
servational signatures, such as type I X-ray bursts or co- 
herent X-ray pulsations, black holes at best offer negative 
evidence: the absence of any clear neutron star character- 
istics. Several criteria have been proposed to distinguish 
black holes from neutron stars in x-ray binaries based on 
solely on their x-ray emission, but so far none has been 
found be truly unique to black holes. It has been proposed 



that two-component spectra, with a thermal component 
and an extended steep power-law component, are a sig- 
nature of the presence of a black hole in a binary system 



(Laurent & Titarchuk 1999). Our data clearly show the 
presence of both a soft thermal component and a steep 
power-law component extending to at least 200 keV from 
an x-ray binary containing a neutron star. Again, a crite- 
rion proposed to distinguish black hole versus neutron star 
binaries based on their x-ray emission is found inadequate. 
The spectrum we observe from 4U 0614-1-091 may also 



be compared to the proposal by Zdziarski et al. (1998) 
that very high energy cutoffs, £^cut > 100 keV, and ther- 
mal Comptonization temperatures, kT^ > 50keV are a 
signature for black holes. However, this criterion applies 
only to the hard state (F < 2) emission, while we found 
4U 0614-1-091 in a softer state (T = 2.45). It would be of 
great interest to observe 4U 0614-f091 with BeppoSAX 
while it is in a harder state to see if the spectrum still 
extends to high energies without a cutoff. 

Our best deconvolution of the spectrum requires re- 
flection from an optically-thick disk. Reflection models 
have worked very well in fitting broad bumps, between 
10-60 keV, visible in the spectra of black h ole candi- 
dates ( Done et al. 1992 ; Gierlihski et al. 1997 ; Zdziarski 



et al. 199§| ). For two other x-ray bursters, 4U 1608 



522 and GS 1826-238, observed by Ginga in their hard 
( low) state the refle c tion model gave also good results 
( Yoshida et al. 1993 ; ^trickman et al. 1996 ). The excel- 
lent broad band coverage of BeppoSAX gives strong evi- 
dence for reflection in the spectrum of 4U0614-f091. We 
find that the magnitude of the reflection and the flux of 
the blackbody component are both correlated with the 



photon index. Zdziarski et al. (1999) found a common cor- 
relation between strength of reflection and spectral slope 
in Seyfert AGNs, 2 black hole candidate X-ray binaries, 
and the X-ray bursters 4U 1608-52 and GS 1826-238. The 
correlation we find for 4U 0614+091 is similar in form but 
may be slightly displaced to higher values of the photon in- 
dex. Additional high-quality broad-band spectra obtained 
with BeppoSAX would permit more detailed study of this 
correlation in the future. 
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